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Abstract. Extensions of the fractional Brownian fields are constructed over a complete 
Riemannian manifold. This construction is carried out for the full range of the Hurst 
parameter a 6 (0, 1). In particular, we establish existence, distributional scaling (self- 
similiarity), stationarity of the increments, and almost sure Holder continuity of sample 
paths. Stationary counterparts to these fields are also constructed. 



1. Introduction 

The fractional Brownian motions and their stationary counterparts are the basic exam- 
ples of Gaussian random fields over R and it is natural to ask what are the corresponding 
examples when R is replaced by a manifold. The first to do so was Paul Levy (see |25j), 
who extended the standard Brownian motion on R to the standard Brownian field over R^, 
now called Levy's Brownian motion. Levy then extended this field to the sphere E> d . Since 
then there have been a number of studies aimed at extending both the Brownian motion 
and the fractional Brownian motion to other manifolds. This is a natural step in the theory 
of Gaussian fields in general as one would like to understand how the structure of the index 
set determines the kinds of fields that can be defined over it. The geometric and topological 
structure of Riemannian manifolds make them a convenient and interesting setting for such 
a study. When one extends the fractional Brownian motions from R to R d the resulting 
fields are called fractional Brownian fields (some authors prefer Levy fractional Brownian 
motions) and our purpose in this article is to construct fields over Riemannian manifolds 
that generalize the fractional Brownian fields over M. d . 

Much of the interest in the fractional Brownian fields (/-B/'s) over R d stems from their 
distributional invariance and scaling properties. In particular, if a G (0, 1) denotes the 
Hurst index and the corresponding field is denoted by fBf a , the increments of the fBf a 
are invariant under rotation and translation and the distribution of the fBf a scales by a 
power c a when M. d is dilated by c > 0. Any extension of the fBfs should possess these 
properties and also reflect the geometry of the index set in question (for an introduction 
to Gaussian random fields over manifolds focusing on smooth fields see the recent work of 
Adler and Taylor [2] and the lecture notes from the short course given at the 2012 joint 
meetings of the AMS^) . 
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As mentioned above the first attempt to extend Levy's Brownian motion, fBf*, from 



d(x, y) being the geodesic distance between x and y and o being a fixed origin point on 
the sphere. Further progress in this direction was made in the work of Molchan (see e.g. 
|28j ) and Gangolli (see [16]) where the authors dealt with extensions of Levy's Brownian 
motion to other manifolds including the sphere. 

Most recently Istas in [21] studied fields over certain Riemannian manifolds with covari- 
ance given by 



where d(x, y) is the metric of the manifold and o is a chosen point. In particular Istas 
showed there that (1.1) defines a Gaussian field over compact rank one symmetric spaces 
and hyperbolic space M d if and only if a 6 (0, 1/2]. 

A common feature of the above approaches is that they begin by looking for covariances 
of the form f(x, 6) + f(y, o) — f(x, y) for some symmetric function /; the idea being that 
over M d o = and f(x,y) = \\x — y\\^d. The issue then is to prove that the function so 
definined does, in fact, define a covariance, i.e., one must establish positive definiteness. A 
necessary and sufficient condition for positive definiteness is that / be of negative type, for 
example one can take the above approach on metric spaces (X, d) with metric of negative 
type (e.g. p2, 20] ). In general if d(x, y) is the metric of a Riemannian manifold, establishing 
that d(x,y) 2a is of negative type for some a G (0, 1) is non-trivial and indeed, as in [21], it 
has been shown d(x, y) 2a can fail to be of negative type. Moreover, in all the above work 
this approach necessitates symmetry assumptions on the underlying manifold. 

In the present article we take an essentially different approach inspired by the work 
of Benassi, Jaffard, and Roux (see [I] and more recently [5]). In particular we extend a 
characterization of the fBf a in terms of the Laplacian on M d to the Riemannian setting 
via the Laplace-Beltrami operator and the associated heat kernel. Using this approach 
we are able to extend the fBf a to a variety of both compact and non-compact manifolds 
without any assumptions regarding symmetry of the manifolds and for the full range of 
a £ (0, 1) (see Theorems 3.2-3.4 below). 

This article is structured as follows: in Section 2 we cover some preliminaries regarding 
Gaussian random fields and analysis on manifolds, in particular the heat kernel of a Rie- 
mannian manifold. In Section 3.1 we describe the motivation behind our approach before 
moving on to discuss existence in Section 3.2. In Sections 3.3 and 3.4 we establish distri- 
butional properties and sample path regularity. In Section 4 we construct stationary coun- 
terparts the fields of Section 3 and establish the corresponding distributional properties. 
Section 5 contains some open questions concerning geometry and probability encountered 
in the course of the article and in the appendix we collect some auxiliary results concerning 
sample path regularity of Gaussian fields over manifolds. 




(1.1) 
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2. Preliminaries 

2.1. Gaussian Random Fields. Given a complete probability space (fi, J 7 , P) and some 
index set / we call a collection of random variables on f2, {Xj(u;)}j g /, a Gaussian random 
field (GRF) over I if for any finite subset {ik}\ C / the random vector (Xj fe )™ has a joint 
normal distribution. Then for each ui 6 f2, -Xj(w) defines a real valued function on I called 
a sample path of the field {Xi}. We let E denote the expectation operator, 



E[Xi]= [ X t (Lu)dP(u) 
Jn 



i G / 

Jn 

and we call 

E[(X S - E[X s ])(X t - E[X t }) = E[X s X t ] - E[X s ]E[X t ] s,t € I 

the covariance of {Xi}. The covariance of a GRF over / defines a symmetric positive 
definite function on I x I. 

We say two GRF's are equal in finite dimensional distribution or simply in distribution, 

denoted =, if their covariances are equal. We also say two GRF's are versions of each 
other if they are equal in distribution. The salient analytical feature of GRF's is that for 
any set / the collection of all GRF's over / is in one to one correspondence up to equality 
in distribution with the set of all symmetric, positive definite functions on I x I. In other 
words a GRF is uniquely determined in distribution by its covariance and every symmetric 
positive definite function on / x / is the covariance of a GRF over /. 

We call a GRF centered if E[Xj] = \/i £ I and in this case its covariance is given by 
E[-Xt-Xs], s,t £ I. Throughout this article we will only consider centered GRF's. 

2.1.1. The Fractional Brownian Fields OverM. d . The standard Brownian motion B t over 
[0, oo ) is the centered GRF with covariance 

^ ^ i \s\ + \t\ — \t — s\ 
E[B s B t ] = s At = U ^A_J -. 

From this one generalizes to obtain the fractional Brownian motion fBm a for a G (0, 1): 

I c 1 2a 4- I i 1 2a — If— c 1 2a 
E[fBm?fBm?} = + '* ■ 



i 



We then have Bt = fBm?. 

One then further generalizes to M d , obtaining the fBf a as the centered GRF over M. d 
with covariance 

E[fBf«fBf«] = \\x\\^ + \\y\\^ d - \\x - y\\^ d 

(note that some authors include the constant factor 1/2). 

One easily sees that the fBf a is self similar of order a, i.e., if fBf® denotes the field 
rescaled field {fBf" x } xeR d then 

fBf? = c a fBf a Vc>0, 
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and that it has stationary (or homogeneous) increments: 

E[|/B/« - fBf«\ 2 ] = \\x - y\\ 2a = \\i(x) - i(y)\\ 2a = E[\fBf^ x) - fBf? {y) \ 2 } 

for any isometry t on M. d . Moreover it is known that there exists a version X x of the fBf a 
such that with probability one the sample paths X x (oS) are Holder continuous of any order 
7 < a and fail to be Holder continuous of any order 7 > a at every point in R d (see P). 

2.1.2. White Noise. The treatment here follows [23J. Given a probability space (Q,J-,P) 
we call a complete subspace G of L 2 (fi, J 7 , P) a Gaussian Hilbert space if every element of 
G is a centered Gaussian random variable. Note that the inner product H inherits from 
L 2 (n,T,P) is then 

(X,Y) G = E[XY]. 

Given any (real) Hilbert space H there exists a Gaussian Hilbert space G and a unitary 
map W : H — > G called W the isonormal process or iw/uie noise process on ff (one can 
also consider complex white noises, but we will only consider the real case). If, as is the 
case below, H = L 2 (M,S,d/j) for some measure space (M,S,dfi) then if B = {A G S : 
n(A) < 00} the map from B — >• G given by 

W(A) = W(xa) 

determines a Gaussian random measure on M. The properties of such measures will not 
be important for us here, but we mention them to motivate the notation for W : H — > G, 
given by 

W(f)= [ f(z)dW(z), 

which we refer to as a white noise integral (this is also commonly called a stochastic 

integral). Starting from a random measure one can construct the integral / dW in close 

Jm 

analogy with classical measure theory. All that will be important for us is the property 

(f,g) H = E [ fdW, [ gdW 
Jm Jm 

Now suppose we have a function h(x, z) : M — >■ L 2 (M, dp), x \-t h(x, z) 6 L 2 (M, dfi(z)). 
We can then define a centered GRF Y x over X by 

Y x = f h(x,z)dW(z). 
Jm 

The covariance of Y x is then given by 

E[Y x Yy] = (h(x,z),h(y,z)) L 2 = / h(x, z)h(y, z) dp(z). 

Jm 

Note that the last expression on the right is in fact positive definite and symmetric. In 
this case we call h the integral kernel of Y. 
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2.2. Analysis on Manifolds. In what follows we assume throughout that all Riemannian 
manifolds are complete and of dimension d, with 2 < d < oo. For a manifold M let A 
denote the Laplace-Beltrami operator, or simply the Laplacian for short, on M. In any 
local coordinate system the action of A on C°°(M) is given by 

where (gij) is the matrix of the Riemannian metric in these coordinates, (g ij ) = (g%j)~ , 

and y/g = (det(gij))^ . Because M is complete, A is essentially self adjoint (see e.g. [50] ) 
and so we may consider from now on the unique minimal self-adjoint extension of A, which 
we shall write as A also. Moreover the spectrum of A is contained in (— oo,0] (see e.g. 
|30j). By the spectral theorem we can define the heat semigroup 

poo 



where dE\ is the spectral measure of -A. The action of e tA on L 2 (M,dV g ), where dV g 
denotes the measure derived from the metric g, is given by a kernel H t (x,y): 



(f)(x)= I H t (x,y)f(y)dV g (y). 

M 



Ht(x,y) is called the heat kernel of M. It is known that Ht is strictly positive, symmetric, 
and contained in C°°(M x M x (0, oo)). Moreover we have the semigroup property 



/ H t (x,z)H s (z,y)dV g (z) = H t+S (x,y). 
Jm 



As a consequence H t is positive definite for each t > 0. As its name suggests, H t (x,y) is a 
fundamental solution to the heat equation on M x (0, oo): 

( (jt-A a \H t (x,y) = 

lim / H t (x,y)f(x)dx = f(y) V/ € C (M). 
L 40 J M 

There are various constructions of the heat kernel, that given in [TUJ being most suited 
to our purposes. In particular if we let 



then there is an open neighborhood of the diagonal U C M x M such that on U 
(2.1) =•(*,,,,) 
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where x, y) is symmetric in x and y, $ G C fc ([0, T] x [/) V T > where fc can be chosen 
arbitrarily large (see [7] and [6]), and 

lim $(t,x,y) = 1. 

I— M), x— >j/ 

In other words, for x and y close Ht ~ £ t as i — )• 0. Thus on any manifold heat diffusion 
behaves locally for small times as in Euclidean space. 

If M is compact then we also have the following eigenfunction expansion of Hf 

oo 

(2.2) H t {x,y) = Y j e- Xkt ^k{x)<P k {y) 

k=0 

where = Ao < Ai < ... < A& f oo and {4>k} are the spectrum and orthonormalized L 2 
eigenfunctions of —A respectively and where (2.2) converges absolutely and uniformly for 
each t > (see [7]). 

Following [7j we define a regular domain to be an open, connected, relatively compact 
subset D of a complete Riemannian manifold such that dD ^ is smooth. In what follows 
when we refer to the Laplacian of a regular domain we mean the Dirichlet Laplacian with 
corresponding the heat kernel (see [7J, Chapter 7). As in the compact case we have an 
eigenfunction expansion (2.2), the only difference being that Ao > 0. If (M,g) is a regular 
domain in manifold (N,g) then, as noted in [8], (2.1) holds in this setting as well. 

Now suppose M is complete and non-compact, {D^}^ is any increasing exhaustion of 
M by regular domains, and H^(x,y) denotes the Dirichlet heat kernel of D^. Then if 
we extend each H k to be zero outside D x D, {H t (x, forms a pointwise increasing 
sequence on M x M x (0, oo). It was shown in [14] that 

lim H?(x,y) = H t (x,y) 

k— >oo 

where H t (x,y) is the heat kernel defined above. 

3. The Riesz Fields 

3.1. Motivation and Definition. In [4] the authors begin by defining a symbol class of 
pseudodifferential operators over M. d . From such an operator A they define a Gaussian 
random field with covariance given by the integral kernel of A^ 1 . The authors are then 
able to derive all the important properties of this field from properties of the symbol of the 
operator A. This approach to constructing and studying GRF's is a natural extension of the 
classical spectral theory of Gaussian processes on R and provides a beautiful demonstration 
of the power of the spectral point of view. 

The basic heuristic can be described as follows: Beginning with an unbounded operator 
A on some L 2 space, define and study the GRF determined by the integral kernel of A~~ l . 
So in attempting to extend the fBf a to a Riemannian manifold, we should first seek an 
operator A that determines the fBf a in the manner above. 
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Our starting point is the well known (e.g. [I] or [32]) spectral representation of the fBf a , 
(3-1) fBft ± C d>a / —r— dW(0, 

where W is a complex white noise on L 2 (M d ,dx), dx is Lebesgue measure, and Cd a is a 

(!+f) ( 

f(0 di = (-A)-(!+f - (-A)-d+f )(/)(0). 



constant. Examining (3.1) we see that, up to a constant, for / £ i2_(d , s \ (M rf ) 



Thus if we denote this last operator above by A then the fBf a is the unique (in distribu- 
tion) GRF with covariance given by the Schwarz kernel of the operator A* A, 

1 yJ J* U\\ d+2a 

We now have a characterization that extends immediately to manifolds: Simply re- 
place the Laplacian on M. d by the Laplace-Beltrami operator of the manifold in question 
and determine the kernel of the operator A* A. Following [30] we arrive at the following 
definitions: 

Definition 3.1. For a complete Riemannian manifold M with heat kernel Ht(x,y) define 
the Riesz field R a to be the GRF with covariance given by 

1 f°° d 

(3.2) nKRy\ = — v / ^ +a ' x ( H *( x > v) - H t( x > °) - H t(y> °) + H ^ °)) dt 

i (2 + a ) ^0 

where o 6 M is a fixed "origin" and the stationary (or homogeneous) Riesz field hR a the 
GRF with covariance 

1 f °° d 

(3.3) E[hR"hRy] = / t*+ a - l H t {x, y) dt. 

1 (2 + a) Jo 

Because H t (x,y) is positive definite for each t > and 
H t (x, y) - H t (x, o) - H t (y, o) + H t (o, o) 

[H t/2 (x, z) - H t/2 (o, z)) (H t/2 (y, z) - H t/2 (o, z)) dV g (z), 



IM 

each of these expressions is symmetric and positive definite, and thus when the integrals 
exist each determines a GRF over M . Of course the convergence of the above integrals is 
by no means obvious and our first task in Section 3.2 will be to determine manifolds for 
which they do converge. 

Remark 3.1. We will see shortly that if either (3.2) or (3.3) exist for some a$ G (0, 1) then 
it also exists for any a 6 (0,ao). We say R a (resp. hR a ) exists for all a G (0,6) if (3.2) 
(resp. (3.3)) is finite for all a E (0, b), b < 1, and all x, y G M. 
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It turns out (Proposition 3.5) that the Riesz field (3.2) extends the fBf a and that they 
agree up to a constant in distribution over M. d . However we will also see that the stationary 
Riesz field has some claim to be an extension of the fBf a , for example over negatively 
curved manifolds, even though it does not exist on M. d . 

Now let W denote the white noise over L 2 (M, dV g ). We will show that when they exist 
the Riesz fields admit the following integral representations: 

(3.4) Ag= / rti^- 1 (H t (x,z)-H t (o,z))dtdW(z) 

1 (4 + 2) Jm Jo 

and 



(3.5) hRzl—Lf r t i^- 1 H t (x,z)dtdW(z). 

1 (4+27 JM JO 

The issue is whether or not the functions appearing in the above are in fact square 
integrable for each x 6 M. Let us consider this in detail, first for hR a : 

1 f°° d 

Letting hhn(x, z) = — j-, / 1 4 + "2 ~ 1 H t (x, z) dt we have 

r (f + f ) Jo 

(h hR (x,z),h hR (y,z)) L 2 

= r ( 1 r t ^- i H t (x,z)dt) 

Jm\T(I + ^Jo J 

x ^^J\^H t (y,z)ds^j dV g (z) 

= I I" I" —r — 2 ti+% ~ lsi+f ~ lH t&z) H s(y,z)dtdsdv g ( z ) 

JmJo Jo r(f + f) 

j-CO />CO 1 d d f 

= / —T* t 4+f-i S 4+f-i H t (x,z)H s (y,z)dV g (z)dtds 

Jo Jo r(f + §) Jm 

= / / —n - 2 -t* + %- 1 s* + %- 1 H t+s (x,y)dtds 

Jo Jo r(i + qy 

poo poo -1 

= / / *(t - s)4 + 2- 1 S 4+2~ 1 ^( x : y)dtds 

Jo Js T(i + qf 
= / / —H, 3 * " s)4 + 2-^-4+ 2- 1 dsH t (x,y)dt 

Jo Jo r(f + f) 

by the positivity of H t (x,y) and the semigroup property. 
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Next note that if q(s) = — -S4+2 1 then 

r(f + f) 

/ At - 5)4^2 A S 4^2 1 ds = g* git) 

Jo T(i + ^f 

where * denotes the finite convolution / * g{t) = I f(t — s)g(s)ds. If C denotes the 

Jo 

Laplace transform we have the well known property £(f * g) = C{f)C{g). Applying this 
to g * g above we have 

C(g * g)(s) = {C{g)f (s) = (*-(2+f ))* = s~^ +a \ 
Then inverting C we obtain 

Substituting this into the integral defining {hf l n(x,z),hf l n(y,z))i2 above yields 



— ^ / t2+ a - l H t {x,y)dt. 

r(| + a)J v ; 



~2 

Thus whenever hR a exists it is given by (3.5). 



1 f°° d 

Turning now to (3.2), let h R (x,z) = — — / t4+f - 1 (H t (x, z) - H t (o, z)) dt. 

r U + 2) Jo 



Then 

|2 



\\h R (x,z)\\i 2 

< s-^-H-^- 1 \H t {x,z)-H t (o,z)\\H s {x,z)-H s {o,z)\dsdtdV g (z) 
jm jo jo 

roo roo r 

= / si+f-44+f- 1 / \H t (x,z)-H t (o,z)\\H s (x,z)-H s (o,z)\dV g (z)dsdt 
Jo Jo Jm 

noo 
g |+f-i t |+f-i|| jfft(a . j _ ^ .)|| 2 || Fs ( X) .) _ Hs ^ .)|| 2 dscft 

= tS+f- 1 ^^, •) - tf t (o, ->ll=2 dt) 

= (J ti+f-^Htix, x) - 2H t (x, o) + H t (o, o) dt) . 

Recall that if M is any Riemannian manifold then from (2.1) for any x,y £ M we have 
that H t (x, y) = 0(t~i) as t ->■ 0. So then 



JO 



t(x, x) — 2H t (x, o) + -£/t(o, o) dt < co 
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and 

•1 



/ 

J 



/o 

for all a G (0, 1). 

Next notice that if a + e < 6 

t4+f-V^t(x,x) - 2ff t (x,o) + H t (p,o)dt 

/oo 

< Q^°° rft) 2 (i^(x, x) - 2H t (x, o) + H t (o, o)) dt^j 2 

by Cauchy-Schwarz. Thus if R a exists for all a G (0, 6) we may interchange the order of 
integration as with hR a to obtain 

(h R (x,z),h R (y,z)) L 2 

1 f'°° d 

= -j- d r / t* +a - 1 (H t (x,y)-H t (x,o)-H t (y,o) + H t (o,o))dt 

r (§ +a) Jo 

= nKRy\ 

for all such a. 

In either case of (3.2) or (3.3) we see that the integrands are continuous on (0, oo) so 
by (2.1) convergence depends only on the behavior of the integrand at infinity. Thus the 
existence of both and hR^. will depend on the large-time asymptotics of Ht(x, y). These 
depend on the manifold in question and we will treat distinct cases below. 

3.2. Existence. 

3.2.1. The Compact Case. We have the following: 

Theorem 3.1. If M is a compact Riemannian manifold, then the Riesz field of order a 
exists over M for any a G (0, 1) and the stationary Riesz field does not exist over M for 
any a G (0, 1). 

Proof. Recall (2.2): 

oo 

H t (x, y ) = J2 e ~ Xkt Mx)Mv)- 

We have 



H t (x,x) - 2H t (o,x) + H t (o,o) = e' Xkt \<p k (^) ~ <Pk{o)\ 2 = 0(e~ Xlt ) Vx G M 

k=i 

and Ai > 0. Then (3.2) is clearly finite for any x G M and all a G (0, 1). 
To see th 

V x, y G M. 



To see that hR® does not exist on M notice that \imH t (x, y) = Vol(M) 1 / 
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□ 

Theorem 3.2. If M is regular domain then hR a , and thus by linearity R a , exists for any 
a £ (0,1). 

Proof. As above let 

oo 

H t (x, y ) = J2 e ~ Xkt Mx)Mv)- 

k=0 

Then A > and H t (x,y) = 0(e~ A °*) for each x,y G M. 

□ 

We note here that in either case above we may integrate term by term using the eigen- 
function expansions of H t to obtain a series expression for the covariance of R a and hR a 
as follows: 

For R a and M compact we have 

E[R«R«] = — / t2 +a - l H t {x, y) - H t (x, o) - H t (y, o) + H t (o, o) dt 

r (I + a) Jo 

1 f°° 00 

r (I + a) J ^ 

r (I + a) 7 ^ 

r (l+«)v^o ~ / 
x (/ t^E^'i^^-^wi 2 *) 

= (E[K| 2 ]E[|i?«| 2 ])^ 

which we know from above to be finite. 

Then by dominated convergence we may integrate term by term to obtain 

1 ~ r(^ + a) 

E[R*Ry] = v(d . x E *+« (<fe(x) - ^(°))^^) - ^*(°)) 
1 I 2 + Q J fc=l A| 

oo 

= E(A fc )-(' +a )w fc (x) - Mo))(Mv) - Mo))- 
k=i 
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The same equality holds for M a regular domain if we number the spectrum as {Afc}^°. 
Similar arguments show that for M a regular domain 



E[hR5hI%] = Y^{\k)-^ +a hk{x)ct> k {y). 



k=l 

Example 3.1. Let M = § 2 . Then in terms of the spherical harmonics {Ifcm} we have 

oo k 

H t (x,y) = Y,e- k{k+1)t Yl Yk m {x)Y km {y). 

k=0 m=—k 

Applying the harmonic addition formula we have 



47T 
k=0 



where P k is the k-th Legendre Polynomial and (x,y) = cos6 xy . Fixing an origin point 
o G § 2 we then have 

oo , 
nKRy] = ^ k + (Pk(cos9 xy ) - Pk(cos9 xo ) - P k (cos6 yo ) + Pfc(l)) . 

fc=l 

Example 3.2. Let M = D = {x G R 2 : \x\ < 1} and J k the Bessel function of the first 
kind of order k, k = 0, 1, 2... Then if < A| < ... are the positive zeroes of Jfc, using polar 
coordinates on D we have 

E^^i?^] = — ^ Vr ahi ^(Air)J fc (Ai^) (cos(fc(0 - </>)) + sm(k(6 + 0))) . 

3.2.2. The Non-Compact Case. For the case of M non-compact, first let us show by ex- 
ample that we cannot establish existence in general. 

Example 3.3. Let M = S 1 x R. Then we have 

H t M ({8, x), (d>, y)) = Hf(6, </>)Hf(x, y) 

where H M is the heat kernel of M, H s is the heat kernel of S 1 , and is the usual heat 
kernel on R (see |19j . Theorem 9.11). 
We then have that 

H^((d,x),(d,x)) -2^ A/ ((^^),(0,y) + H t M ((M,(M) ' 



7T ^/(Srt) 
3 

= O(ta) as i — >• oo 

for any (0,x), ((f), y) G Af. So E^] 2 ] = oo Vp G M and a > 1/2 and thus R a does not 
exist over M for this range of a. Using § 2 instead in the above we have that R a fails to 
exist for all a G (0, 1). 
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However, for certain manifolds such that Vol(M) < oo we have a situation similar to the 
compact case: 

Theorem 3.3. Suppose M is non-compact with Ric(M) > —k 2 , k£|, and Vol(M) < oo. 
Let A(M) = inf {Ai : cr(S7) = {A^I^q} where the infimum is taken over regular domains 

f2 C M and <r(f2) denotes the Dirichlet spectrum ofO,. Then if A(M) > R a exists over 
M for any a G (0, 1) and hR a does not. 

Proof. That hR a does not exist follows from the fact that on such M 
For R a , under the hypothesis of the theorem it was shown in [23] that 

1 / \(M) 

Htix > y) -voW) =0 { e ~ 2 ' 

and so (3.2) converges V a G (0, 1). 

□ 

We now turn to our main existence theorem for the Riesz fields over non-compact man- 
ifolds followed by some examples. Below we use the following notation: 

D p( r ) = {x e M : d(x,p) < r} 

and 

V p (r)= Vol (D p (r))= [ dV g . 

JD p {r) 

We write Ht = 0(t~z) if there exist two distinct points Xk 6 M, k = 1,2, and constants 
Ck > such that 

H t (x k ,Xk) <C7 fc t-3 Vt>l. 
In that case using Theorem 1.1 of [18] we know that for any 5 > there exists a constant 
Cs > such that for all t > 1 and all x,y E M 



H t (x,y) < C s t~ie ( 4 +*)* . 

Theorem 3.4. Let M be non-compact. 

(1) Suppose Ric(M) > 0. Then hR a does not exist for any a G (0, 1). If 



and 

y ( r ) 

lim -p^ < oo VxeM 

r— >oo r " 

for some (3 G [0, 1) then R a exists over M for any a G (0, 1 — f3). 
(2) Suppose that 



for some /3 > 0. Then hR a (and thus R a also) exists for any a G (0, min{/3, 1}). 
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-(d/2-/J) 



Proof. (1): To begin we note that our hypothesis Ht = 0(t ^ ' p ') implies the following 
gradient bound for Ht (see [H]): For all x,y £ M and t > 1 



(3.6) \V x H t (x,y)\ <C' 5 t~d-^)e 



1 \ d(x,y) 
(4+<5)t 



for some constant C' s > 0. 

Recall that in order for for (3.2) to converge it is sufficient to show that 



oo 

t % +o-l (H t (x, x) - 2H t (x, o) + H t (o, o)) dt < oo 

for the specified range of a. To that end let D = D p (r) be a disk containing x and o. We 
first apply the mean value theorem: 



/oo 
1 (H t (x, x) - 2H t (x, o) + H t (o, o)) dt 



00 d 



Vr**- 1 / |# t 0r,z)-H t ( O ,*)|W s (z)di 

<d(x,o) 2 r^- 1 [ \V x H t (£ z ,z)\ 2 dV g (z)dt 
Ji Jm 

for some £ 2 lying on some curve (parametrized to have unit velocity) contained in D p and 
joining x and o. We now apply (3.6), 



We have 



Ji Jm 



fOO 

<C \ t -2 +a+2 ^ 2 / e ~TW" dV g {z)dt. 

M 



t - 1+ a+2P-2 I e -^+Sjt dVg(z) dt 

Jm 

d +a+2p-2 



t-^+a+ip-i / e (4+4)t dVg ^ dt 
JD 

+ r t -S+o+V>-2 f e - 2J ^dV g (z)dt 
Ji Jm\d 
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/oo d 
t~2 +a+2 P~ 2 dt < oo so we only need to show 

r r 00 2d(£ z z) 2 

/ / t-z +a+2p - 2 e--&riv r dtdV g (z) < oo. 
Jm\d Jo 



We have 

f-OC 



r r t -i +a+ w-2 e -*-f^ dtdVg{z) 

Jm\d Jo 

= ( i ^) 5 " a " 2/3+lr (^-«- 2 / 3 + 1 ) f M ^d(^zT d+2a+ ^- 2 dV g {z). 

Recall D = D p (r) and let 

Ak = D p (r + k)\D p (r + k — 1) k = 1,2,3... 
By monotone convergence 

/ dfe, g )-»**W-* dV g {z) = V / *)-* f2o4 ^- a 

JM\fl A . =1 J A k 

Vol(A k ) 



— ( r _|_ _ l^d-2a-4P+2 

= A ^ P (r + fe) - ^ P (r + fe - 1) 

~~ ( r _|_ ^ _ ]_)d-2a-4/3+2 



Because Ric(M) > we have (see [12] or [9]) 

Vp(cr) < c d V p {r) Vr>0, c > 1. 

Thus 

F p (r + fe) - y p (r + fe - 1) ^ M r + k-l) 



(r+k) d -(r+k-l) d 

EV p yi T K) — V p (l T K — s.) \ - i y (r+k-l) d 

( r _|_ _ l^d-2a-4/3+2 — ( r _|_ J. _ ;[)d-2a-4/3+2 

(r+fe-l) d 



fc=l 



( r + fe - l)^- 2 / 3 f i <■+/■■ r-., -/,-[,■■ 



OO 



(r + k) d - (r + fe - 1) 

Z-/ ( r _|_ £ _ ]^d-2a-2/3+2 



d 



The convergence of this last sum is equivalent to that of 

00 h.d-1 00 

E K = V jt.2a+2/3-3 

£,c2-2a-2/3+2 Z-/ 
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By hypothesis a < 1 — ft, which implies 

oo oo 
k=l k=l 

for some e > 0. 

To see that hR a does not exist on M for any a, we note that (see e.g. |12j ) 

Ric(M) > H t (x,y) > (4vrt)~i e'^r" 
for all x,y £ M and t > 0. Thus 

rco 



/•oo 

/ t2 +a - l H t {x,y)dt = oo 
Jo 



for all x,y £ M and any a G (0, 1). 
To prove (2), simply write 



/oo roo 
1 2+ a ~ 1 F i (x, y)dt<C J t a -P~ l dt < oo. 



□ 



Remark 3.2. There is by now a large body of literature relating geometric or functional 
analytic conditions on M and the large-time decay of the heat kernel. We will not attempt 
to make a survey, but only mention the names of some researchers that have worked in the 
area recently: I. Chavel and E.A. Feldman, T. Coulhon, A. Grigor'yan, and E.B. Davies, 
among others. We suggest that the interested reader begin with the work of these authors 
and the references therein. 

We are now in a position to show that, over M. d , R a agrees up to a constant with the 
fBf a in distribution. We could do this abstractly using arguments along the lines of 
Section 3.1, however we can also make a simple explicit calculation. Note that M d satisfies 
the first hypothesis of Theorem 3.4 with (3 = 0. Thus R a exists there and if we choose 
o = has covariance 

E{R«R°] = — d / t2 +a -\H t (0,0) - H t (x,0) - H t (y,0) + H t (x,y))dt. 

, 1 -|k-y|l 2 j 

Proposition 3.5. If M = R then HAx.y) = — e « and for all and 

for a £ (0, 1) 

E[RZR%] = C a (\\x\\ 2a + \\y\\ 2a - \\x - y\\ 2a ) 
where C a is the positive constant given by 

-r(-a) 



C' 



42+ a (vr)2r(f + q) 
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Proof. First note that if either x = or y = the result is trivial; thus we assume otherwise. 
The integral defining E[R%Ry] reduces to 

°° „ i/-, Hl^ll 2 -Nil 2 -lk-»ll 2 , , 



A/(47r) rf Jo 



— M J 'II — 1 1 y 1 1 — II — 

which we recognize as a Mellin transform. Let F a (t) = X[a,oo){t) ~ e 4t — e « + e «~ 
with a > 0. Then F a {t) = 0{t) as t oo and F a (i) = o(t Ar ) as t V iV > 0. Thus 



poo 

/ t s - l F a {t)dt 
Jo 



10 

converges absolutely for all s G C with !R(s) < 1 and defines an analytic function there. 
On the other hand for — 1 < K(s) < we have by direct calculation that 

f 00 , lnM1 a s -||x|| 2s - ||y|| 2s + \\x- y\\ 2 ° . . 

By analytic continuation this last equality holds for < 3?(s) < 1 as well. For such s we 
have by dominated convergence 

f t s - l F Q (t)dt = lim f t^Fai^dt. 
Jo Jo 

Now for a < 1 

/OO /'OO 
t s - 1 F a (t)dt = J t s - l F Q (t)dt 

and so, noting Fq{€) > 0, we have using dominated convergence 

poo /"I poo 

/ t s - 1 F (t)dt= / t s - 1 F (t)dt+ / t s ~ l F (t)dt 
Jo Jo Jl 

= ( lim / t^Fa^dt] + f t'- 1 F (t)dt 

\a^0+ J J Ji 

= lim f / t s ~ 1 F a (t)dt + l°° t^F^t) 

a^0+ \Jo Jl 

pOO 

= lim / t s ^F a (t)dt 

a^0+ Jo 

= ~IMI 2s -\\y\\ 2s + \\x -y\\ 2s F ( s j 



□ 



Example 3.4. Suppose M is non-compact with Ric(M) > and 

lim ^ = ^(0,1) 
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for some p E M (cf. the Bishop-Gromov comparison theorem). Then R a exists over M for 

any a £ (0, 1) and hR a does not. Indeed, in [26] it is shown that H t (x,y) = 0{t^^) for 
every x,y 6 M. Theorem 3.4 applies once we note that for all p G M 

Ric(M) > V p (i?) < V > 0, 

ujd being the volume of the unit ball in M. d . 

Example 3.5. If M is simply connected with all sectional curvatures K < k for some 
k < and Ric(M) > —k 2 > — oo then exists over M for any a > 0. For example 
this holds if M = M d , ci-dimensional hyperbolic space. This follows from [27] in which it is 

\k\ 

shown that a (—A)) C [(n — 1) — , oo), which in turn implies the following upper bound 
on H t (see [13]): 

(n-l) 2 fct 

H t (x,y)<Ce — 2 — Vt > 1 

for some C > and all x,y £ M. Theorem 3.4 then applies. 

In particular for M = H 2 , letting p = d(x, y) we have the well known formula 

H t (x,y) = - -3-e * / — — w-rds- 

(4vrt)2 J P cosh(s) - cosh(p) 

Then 

E[fci^fci2?] = 3^ / / i a "3 —dsdt. 

(4vr) 2T (1 + a) cosh(s) - cosh(p) 

Remark 3.3. On negatively curved manifolds, can also be viewed as an extension of 
the fBf a in the following way: In Section 3.1 we saw how the covariance of the fBf a is 
the integral kernel of the operator A* A where 

A(f) = (-A)-(2+f>(/)(:r) - (-A)-(l+f)(/)(0), 

which can be seen as a correction to (— A^^+t) when this operator does not have an 
integral kernel. However on manifolds with spectrum as in example 3.5 (— A) _ *-4+2') does 
have an integral kernel and no correction is needed. So if we view the fBf a as the GRF 
with covariance that is the kernel of the minimal correction to (— A)~^4 + 2) that yields an 
integral operator, then on such manifolds as above we obtain the hR a . 

3.3. Distributional Scaling and Invariance. Having established existence on a variety 
of manifolds we turn next to establishing distributional and sample path properties of the 
Riesz fields. 

Definition 3.2. Let (M,g) be a complete Riemannian manifold and I{M) the group of 
isometries of (M,g). If Y x is a centered GRF over (M,g) we say that Y x is stationary (or 
homogeneous) if 

E[Y l{x) Y L{y) } = E[Y x Y y ] 
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for any t G I(M) and all x,y £ M. We say Y x has stationary (or homogeneous) increments 
if 

ny,{ X )-y,(y)\ 2 ] = n\y,-y y \ 2 ] 

for any i G I(M) and all x,y £ M. 

Because for any manifold (M,g) we have Ht(c(x), i(y)) = Ht(x,y) for any l £ I(M) (see 
[19J, Theorem 9.12) it is clear from the definitions, 

1 [°° d 

nKRy\ = r / t^- 1 (H t (x, y) - H t (x, o) - H t (y, 6) + H t (o, o)) dt 

r (§ + a) Jo 

and 

1 f°° d 

nhR a x hR a y ] = -r d ; r / t2+^H t (x, y) dt, 

that when they exist, i? a and hR a have stationary increments and are stationary respec- 
tively. 

If the Riesz fields are to generalize the fractional Brownian fields they should possess a 
distributional scaling property analogous to that of the fractional Brownian fields over R^. 
In order to show this we must first define what it means to scale our manifold (M, g) by a 
factor c > 0. By scaling we mean changing the metric g on M so that d(x, y) (— )■ cd(x, y) 
and the way to do this is to replace g by c 2 g. If (M, g) is a regular domain in (N, g) then 
to scale M we scale iV and so the following applies as well to regular domains. 

Proposition 3.6. Let (M,g) be a complete Riemannian manifold or regular domain. Both 
the Riesz field R a and the stationary Riesz field hR a over (M, g) are self-similar of order 
a (if they exist on M ) in the sense that if R a and hR a are the Riesz fields over (M, c 2 g) 
then 

and 

c a hR x = hR x 

for any c > 0. 

This clearly agrees with the usual definition of self-similarity for a random field Y x over 
R d , c a Y lx = Y x . 

c 

Proof. First we note from the coordinate expression for A, if we denote by A g the Laplacian 
of (M,g) and H^(x,y) the corresponding heat kernel, we have A c 2 g = — A g . But then 



20 ZACHARY A. GELBAUM 

because L 2 (M, dV g ) = L 2 (M,dV c 2g) we can write 

/ c d Hf g (x,y)f(y)dV g (y)= [ Hf g (x,y)f(y)dV c i g (y) 

JM JM 

= e~ tA ^(f) 
= e-£ A »(/) 

= / H 9 t (x,y)f(y)dV g (y) 

JM C* 

for any / G L 2 (M, dV g ). Thus by symmetry 

1 2 

-jH g t (x,y) = Ht 9 (x,y) Vi,y£M. 

We then have 

„2a /-oo 

y r(§ + a)y 

= F7 ^r^ 1 (flf'(x,y)-flf"( 0) x)- J Hf'( 0| y)+flf*( , )) A 

r(§ + a)y v / 

and similarly for hR a . 

□ 

Remark 3.4. Here we see that hR a exhibits essentially non-Euclidean phenomena; on M. d 
there cannot exist a GRF that is both stationary and self similar (see e.g. [3]). We will 
return to the questions this raises in Section 5. 

3.4. Holder Regularity. If M is any Riemannian manifold or regular domain with heat 
kernel Ht(x,y) then the maximum principle implies 

H t (x,y) < H t (x,x) Vx,yGM 

with equality if and only if y = x. We then have that 

H t (x,x)-2H t (x,y) + H t (y,y) >0 Vy^x. 

In particular E[|i?° — R g \ 2 ] and E[|/ii2" — hRy\ 2 ] both define metrics on M when they exist. 
Note also that 

nK-Ry\ 2 \=nhK-hRy*\ 2 \ 

when both exist. In particular in the proof below we will not distinguish these two metrics 
as the context of the Theorem will make clear which is being discussed. 
We are now in a position to prove the following: 



/ 

Jo 
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Theorem 3.7. Let M be a compact Riemannian manifold, a regular domain, or non- 
compact under the hypothesis of Theorem 3.4- Then R a (resp. hR a ) has a version such 
that with probability 1 all sample paths are Holder continuous of order 7 for any 7 < a and 
fail to satisfy a Holder condition on a dense subset of M for any 7 > a. 

Proof. In order to apply Theorem A. 4 in the appendix we need to compare the metric 
E[|i2£ - Ry | 2 ] (resp. E[\hR% - hR^\ 2 ]) on (M, g) with the metric d(x,y) derived from g, in 
particular study the boundedness of 

nR a x -R a y \ 2 } 

K ' (d(x,y))^ 

for d(x, y) small and 7 6 (0, 1). What we will show is that this ratio is unbounded if 7 > a 
and approaches zero if 7 < a. 

Our approach to controlling (3.7) will be to split the integral defining E[|i?° — Ry\ 2 ] into 
two parts: 

00 

t-2 +a ' l (H t (x, x) - 2H t (x, y) + H t (y, y)) dt 

'0 

(3.8) = / ti +a ~\H t (x,x)-2H t {x,y)+H t (y,y))dt 

Jo 

/oo 
t* +a -\H t {x,x) - 2H t (x,y) + H t (y,y))dt. 

We start with (3.8). Recall that in any case around any point p there is a closed disk 
D p such that (2.1) holds with G C k (D p x D p x [0, T]) where we can choose k > 2 and 
T>0. 

As a consequence we have, denoting the integral (3.8) by I\ and d(x, y) by p, 
(3.10) I 1 = (4 ir )-$ [ t Q - 1 ($(t,x,x) + $(t,y,y) -2$(i,x,y)e-£)dt. 

JO 

Because $ G C k (D p x D p x [0, T]) with k > 2 and is symmetric, by Lemma A.l in the 
appendix, 

$(t, x, x) + $(t, y, y) - 2$(t, x, y) = 0{p 2 ) as p -> 
uniformly for t G [0, 1]. Thus we have 

t a_1 (^(t, + i/, y) " 2$(t, x, y)e- « ) (it = ^-^(l - e~£)dt + 0{p 2 ). 
Jo Jo 

Because 

f t a -\l-e~£)dt = p 2a f r-\l-e-*)dt 
Jo Jo 



(■OO 

and 

'o 

(3.11) h = 0(p 2a ) = 0(d(x, y) 2a ) as d(x, y) -»• 



r 00 i 

/ t a - 1 (l- e -«)(it<oo, 
■/ o 
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for x, y G D p . 

For (3.9), which we denote I2, we first deal with the case of M compact. Using (2.2) we 
have for t > 1 

00 

s) - 2H t (x, y) + F t (y, y) = J] e~ Afc *|^(x) - </> fc (y)| 2 

fc=0 

= 2>- Ah V*(*)-fo(v)l 2 

fc=l 

00 

<ci(x,y) 2 ^e- Afct ||V0 fc || oo . 
fc=l 

Now we apply the following bound on || V0fc||oo (see [21?]): 

||V0 fc ||oo < CmA^ 1 

where Cm is a constant depending only on M. We then have 

00 d 

H t (x, x) - 2H t (x, y) + H t (y, y) < C M d(x, y) 2 ^ e^X^ = C M d(x, yfO , 

fc=l 

which yields 

(3.12) I 2 < C M d(x, y) 2 J™ ti +a - l O (e _Alt ) dt = Cd(x, y) 2 

as Ai > 0. 

If M is a regular domain then a similar argument using the corresponding bound (see 

for the Dirichlet eigenfunctions on M we obtain (3.13) in this case as well. Thus for either 
M compact or a regular domain 

h = (d(x, y) 2 ) as d(x, y) -> 0. 

Turning now to the case of M non-compact, first suppose the first hypothesis of Theorem 
3.4 is in force. As in that proof we have, for x,y contained in a sufficiently small geodesic 
disc, 

/oo 
t 2+"-l ( Ht (x, x) - 2H t (x, y) + H t (y, y)) dt 

= ' t l+ a - 1 f \H t (x,z)-H t (y,z)\ 2 dV g (z)dt 

Jl J M 

<d(x,y) 2 / %f+ Q - 1 f \V x H t (£ z ,z)\ 2 dV g (z)dt, 

Jl J M 

which was shown to be finite. 
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Next suppose the second hypothesis holds. For this case we will use a Schauder esti- 
mate and Lemma A.l: We choose a geodesic disc D p and let L be A in geodesic normal 
coordinates on D p , D = exp~ l (D p ) 1 P = dt — L on C°°(D x (0,1)), and u(x',y',t) G 
C°°(D x D x (0, 1)) be Ht(x, y) in our chosen coordinates. For any T > we then have 

Pu(x', y', t + T) = d t u(x', y', t + T)- L x ,u(x', y',t + T) = 

for each for all x ,y ,t E D x D x (0, 1/2). In other words, u satisfies Pu(x ,y ,t) = on 
D x (T, T + 1/2) for each y' G D and T > 0. 

Because L is uniformly elliptic on D and its coefficients are all C°° (and independent of 
T, t), using the Schauder estimate (Theorem 5 p. 64 in [15] and choosing a = 1) we obtain 
for each closed disk D r contained in D a constant K r > such that 



sup 

O',t)eD r x(0,l/2) 



2 



ii 



— {x,y ,t + T) 



dx^x'j 



< K r sup \u(x , y , t + T)\ 
OV)eD r x(o,i) 



for each i,j and y € D r . We then have 



sup 

(i',!/',l)eD r xD r x (0,1/2) 



d 2 u 



— {x,y ,t + T) 



dx'tx'j 



< K r sup \u(x , y , t + T)\. 

(i',j',t)gD r xD r x (0,1/2) 



We note that K r is independent of T and by our hypothesis 

sup Ht(x,y) < Ct~(z + P\ /3 > 0. Thus, returning to D r = exp(D r ), for all T > 1 

(i,!/)eDpXDp 

(x,y,t + T) 



sup 

(x,j/,t)eD r xD r x(0,l/2) 



Then applying Lemma A.l and assuming without loss of generality we have chosen our 
disc Dp such that the above estimates hold, we have 



i; 



f2 



+«-l 



H t (x, x) - 2H t (x, y) + H t (y, y)) dt 



<Cd(x,y) 2 
<Cd(x,y) 2 



+OL-1 



sup 



E 

d 



d 2 H 

dxidx, 



(t,x,y) 



dt 



ti +a -\t-l/2)-^dt 



t2 +a - l (t - l/2) _( 2+^) dt < oo. Lastly recall 
that when hR a exists for a £ (0, b) for some b > then R a does as well. Moreover in that 



case 



nK-Ry\ 2 \ = n\hR a x -hR a y \\ 

so in the second case of Theorem 3.4 the arguments above apply to R a as well. 
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Thus in each case from the preceeding discussion we know that for each p G M there 
exists a closed disc D p centered at p such that for all 7 < a 

n\K-Ry\ 2 ]<c p {d^y) 21 ) 

for some constant C p > (or similarly for hR a ) and all x, y G D p and that such a condition 
fails for any 7 > a in light of (3.11). Then by Theorem A. 4 in the appendix R a (resp. hR a ) 
is almost surely uniformly Holder continuous over D p of order 7 for any 7 < a. Moreover 
from the discussion following Theorem A. 4 there is a dense subset of D p on which R a (resp. 
hR a ) fails to satisfy any Holder condition of order 7 > a with probability 1. By modifying 
R a (resp. hR a ) on a set of P- measure zero we then obtain the same result for any compact 
K C M. Lastly by taking a compact exhaustion of M we obtain the statement of the 
Theorem. □ 



4. The Bessel Field 

We now turn to constructing stationary counterparts to R a by analogy with the Brow- 
nian motion and Ornstein-Uhlenbeck processes on R. We define the Bessel Field of order 
a G (0, 1) by 

(4.1) B«± \ [ r t ^- 1 e- t H t (x,z)dtdW(z), 

which extends the Ornstien-Uhlenbeck fields with covariance given (up to a constant) by 



J 



- — <ie. 



These fields are altogether more well behaved than the Riesz fields, which is not surprising 
in light of the analogy with the Riesz and Bessel potentials. 

Theorem 4.1. The Bessel field exists over any complete Riemannian manifold or regular 
domain M for all a G (0, 1). 

Proof Proceeding as for hR a , for each x,y G M 

E[ggB»]= ( \ ) [ rti+^ 1 e- t H t (x,z)dt [°° si^- 1 e- s H s (y,z)dsdV g (z) 
\i (4 + 2) J Jm Jo Jo 

-( F (jVi))Yr <<+i - i - <+, - ie " Mff - (i -'' ) *' b 

1 f°° d 

r (f + a) Jo 

From the fact that the heat kernel always satisfies lim H t (x,y) < 00 for any x and y, we 

t— >oo 

see that (4.2) converges everywhere on M x M. 

□ 
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Clearly B% is stationary and we can see that it does not possess the scaling properties 
of the Riesz fields. Turning to sample path regularity we have the following result. 

Theorem 4.2. The Bessel field B a has a version with sample paths almost surely uniformly 
Holder continuous of order 7 for any 7 < a and almost surely failing to satisfy a Holder 
condition of order 7 for any 7 > a on a dense subset of M . 

Proof. Split the integral 

1 f°° d 

E[\B«-B«\ 2 ] = — d / t* +a2 - 1 e- t (H t (x,x)-2H t (x,y) + H t (y,y))dt 

1 (2 + a ) Jo 

: -(/1 + /2) 



r(| + a) 



where 



and 



Jo 



h = I t 5+«2-i e -* (H t (x, x) - 2H t (x, y) + H t (y, y)) dt 



h = J ^+« 2 -! e -* (H t (x, x) - 2H t (x, y) + H t (y, y)) dt 
and argue as in Theorem 3.7. 



□ 



5. Conclusion and Further Work 

Using a spectral theoretic approach we have constructed analogues of the fractional 
Brownian fields over arbitrary compact manifolds and a wide class of non-compact man- 
ifolds. However we did not establish uniqueness and, at least for a G (0,1/2), it seems 
unclear what such a result would look like in light of the fact that there now exist two 
distinct GRF's over S rf and M d that have stationary increments and are self-similar. Thus 
one could ask how many different such fields there are over any given manifold. 

We also saw in example 3.3 that R a does not exist on S 1 x R for a > 1/2. This raises 
the following question: Does there exist any Gaussian field over S 1 x R with stationary 
increments that is also self similar of order a for some a £ (1/2, 1)? More generally, are 
there geometric conditions that ensure a given manifold can have such a field defined over 
it? 

There is one aspect of this theory we did not touch upon, that being the behavior of 
our fields when restricted to geodesies. For example one could require that a field when 
restricted to infinite geodesies became a fractional Brownian motion. However this would 
require an infinite geodesic and a family of isometries of the whole manifold that restrict 
to translation of the given geodesic, which is a strong condition indeed. 

We also mentioned above that the existence of hR a raises interesting questions regarding 
negatively curved manifolds and what we could loosely call hyperbolic Gaussian random 
fields. For example, although the proof of existence of hR a over M d is spectral theoretic, 
one can ask if there are more geometric or topological conditions one can put on a manifold 
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M to ensure the existence of some self-similar and stationary GRF. Conversely one can 
ask what are the implications of such a field existing over M. Is hR a the only such field 
or are there others? 

One can view a GRF over a manifold as a kind of "randomization" of the manifold. One 
can then roughly summarize the above questions in the following way: Given a manifold 
M, in what ways can we randomize M and how does the answer depend on the geometry 



The above is only a first attempt to state some questions at the intersection of geometry 
and probability that, at least on the face of it, seem novel and interesting; doubtless there 
are others. The study of random fields over manifolds, although its history is not short, 
seems to the author to still be wide open. It is our hope that the work here and the 
questions raised above will be of interest to both geometers and probabilists and lead to 
further interaction between the two. 



The author wishes to express deep gratitude to his advisor, Harold Parks, for his confi- 
dence, generosity, and encouragement. The author also thanks his teachers David Finch, 
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First we record the following Lemma involving Taylor approximation. 

Lemma A.l. Let M be complete and suppose f G C°°{M x M) is symmetric. Around 
any point p G M there exists a closed geodesic disk D p centered at p and a constant C p > 
such that 



for all x, y £ D p . 

Proof. Let F £ C 2 (R d ) and recall Taylor's Theorem: for each p £ R d and all x £R' 



of Ml 
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Appendix A. 



d 



d 2 f 



\f(x,x)-2f(x,y) + f(y,y)\<C p d(x,y) 2 sup 

DpXD p i - =1 



dxidxj 
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Now let f eC 



2/W 



and f(x,y) = f(y,x). Fix x,y G M. d . Then letting p = (x,y), 



from the symmetry of / we have 
f(x,x)-2f(x,y) + f{y,y) 



d ,.i 

Y\{xi-yi){xj-yj) I (l-t) 

• • -, JO 

d „i 

+ yZ(xi-yi)(xj -yj) (l-t 
■ i JO 

70 ij=l 



d 2 / 

dxidxj 



(x + i(y — x),x)dt 



d 2 f 
dxidxj 

d 2 f 

dxidxj 



(y + t(x -y),y)dt 



(x + t(y - 



+ 



a 2 / 

dxidxj 



(y + t(x - y),y) dt 



c^2(xi- yi)(xj - yj) 



d 2 f d 2 f 
-K-^-ix + 0i, x) + o-My + ^2,y) 



for some constant c > and G M rf with ||#fc||Rd < \\x — j/[| rc j . In particular for x,y in a 
closed disk iX of radius e > we have 



\f(x,x) - 2f(x,y) + f(y,y)\ < Ci\\x - y\\ Rd sup 

D E xD e 



E 



dxidxj 

1,3=1 J 



for some C\ > 0. 

Now suppose / G C°°(M x M) is symmetric and let -D p be a geodesic disk centered at 
p G M. Then the above implies 



(A.l) 
for all x, y G Z? p . 



- 2f(x,y) + f(y,y)\ < C 2 d(x, y) 2 _sup 

D v xD v 



E 



3 2 / 



dxidxj 



□ 



A.l. Continuity of Gaussian random fields. Here we provide analogues of results 
given for Gaussian fields over M d in the setting of manifolds. These proofs are simple 
modifications of the originals and we include them for convenience. The first result is an 
analytical lemma, given for hypercubes in Mr. We will replace the cubes with metric disks 
and M d by a d-dimensional manifold M. Let p be even and continuous on [—1, 1], p(|x|) 
monotone increasing, and satisfy lim p(x) = 0. 

x—>0 



Lemma A. 2. (Manifold version of Lemma 1 in |17j ): Let f G C(Iq) where Lq C M is 
a closed metric disk Iq = {y G M : d(x,y) < r} contained in some normal neighborhood 
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about x £ M. Suppose that 

/ / ex P\ ~TJ- 7tty\ dxdy<B 

JdJd \p{diam{D))J 

for all closed metric disks D C I® . Then for some C > 

\f(x) - f(y)\ < y/]og{BCu-**)dp{u) 

for all x,y £ Iq. 

Proof. Fix x, y £ Io and let Dq be the smallest metric disk containing both (so that 
diam(Do) = d(x,y)). Then choose a sequence of disks D k that shrink to x and such that 
if d k = diam(Dfc) we have 

P{dk) = -p(d k -i). 



exp 



Let fo k = Vol(Ufc) 1 / f dV. We apply Jensen's inequality to obtain 

fp k ~ fDk-A 2 
p(dk-i) J 

< [VoK^OVolp^)]- 1 / / exp( f{x) ~ f [ y) ) dV(x)dV(y) 



<D k JD k _ 1 

< B[yol(£> fc )Vol(£> fc _i)] _1 - 
We then have 



(A.2) \f Dk - / Dfc _J < p(4-i)Vlog(5[Vol( J D fe )Vol( J D fe _ 1 )]- 1 ) 

By the definition of we have 

p(4-i) = 4[p(<4) -p(4+i)] 

and by hypothesis on Jo, 3 C > such that 

Vol(I> fc ) > C(d k ) d 

so that d fc+ i <u<d k ^> u~ 2d < C[Vol(Z>fc)Vol(.Dfc_i)] -1 . Then we can write (4.1) as 

\f Dk ~ /D fc _J < 4 f k ^log(BCu-™)dp(u). 
Summing these and using continuity of / we get 



1/0*0 - fo \ = hm \ f Dk - f Do \ < 4 / Jlog(BC«- M ) dp(«). 

k^oo J V 

Now di < d(x, y) so if we need to we can replace B by a larger bound to ensure the 
integrand is defined, and after doing so we have 

\f(x) - f Do \ < ±J^ V) y/log(BCu-^)dp(u). 
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The argument is symmetric in x and y, so an application of the triangle inequality yields 
the conclusion. 

□ 

Suppose now we are given a (centered) Gaussian random field X x over (M,g) and 
consider its restriction to a closed disk Iq as above. Suppose further that the function 
K(x,y) = K[X x X y ] is continuous on Iq x Iq. Then K(x,y) determines a positive trace 
class integral operator on L 2 (Io,dV g ) and by Mercer's theorem we have 



K {x,y) = ^2 >^k4>k(x)<Pk(y) 



k=0 



uniformly on Iq x Iq, where and <f>\. are the eigenvalues and eigenfunctions of K respec- 
tively. 



Let p(u) = sup{yK[\X x — X y \ 2 ] : d(x,y) < \u\} and X™ = \f^k~4>k( x )@k where the 

k=0 

9k are independent standard normal random variables. 

We then have the following adaptation of Garsia's theorem to the manifold setting: 

Theorem A. 3. (Manifold version of Theorem 1 in |17| Suppose that for x,y £ Iq as 
above 

diam(l ) 



a/— log(-u) dp(u) < oo. 



o 



Then with probability 1 



d(x,y) 



\X? - X™\ < C \^M(Bjp(d(x,y)) + jf y/-]og(u)dp(u) 
where C > and 

f f 1 (X™ -X™\ 2 

B = sup / / exp- dV(x) dV(y) 

mJj Ji 4 \p(d x,y) J V ' W 



is almost surely finite. In particular the partial sums X™ are almost-surely equicontinuous 
and uniformly convergent on Iq. 

Proof. Let 

^- eXP 8{p(d(x,y))J -^ ex P 8 { p (d(x,y)) J 
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where Y k (x,y) = \f\k(<j>k{x) — 4>k{y))6k- Then by independence of the 6k and Jensen's 
inequality for conditional expectation 

E[-P n +1 \Pn, Pi] 

■ -rp- I r ; - I Pn, ■■■■ P\ 



= Pn I E 



> P n exp^ ( E 
= P n exp- I E 



8 ^ p(d(x,y)) 

(Y n+1 (x,y)) 2 - 2Y n+l (x,y)(X™- 1 - X n ~ 1 {y)) 



(Y n +\x,y)Y 
p(d(x,y)) 



p(d(x,y)) 



Pn, Pi 



P 



ni ■ • ■ 1 * 1 



>Pn 



a.s. 



Thus {P n } is a submartingale. Next note that E[P^] < \/2, as 

yn yn 
A i — A y 

p(d(x,y)) 

is centered, Gaussian, and has variance less than or equal to one. Then applying the 
classical submartingale inequalities we have 

EfmaxP^J < 4E[P^] < 4>/2. 

m<n 

Applying the Fubini-Tonelli theorem we then have 



E 



i fxz-x. 



, . max (./•/)- . 

JloJlo m<n 4 \p(d{x,y)) i 

Letting n tend to infinity and applying monotone converge yields 

E[B] <4V2(V(I )) 2 < oo. 
We then have that almost surely 



n \ 2 



£ dV(x)dV(y) <AV2{V{I Q )f. 



so that Lemma A. 2 applies. 
Lastly note that from 



ifxi_ 



X. 



n \ 2 



dV(x) dV{y) < B < oo Vn 



E 



^2 Xk0 k 



.k=0 



oo „ 

^A fc = / K(x,x) dV(x) < oo 
fc=o 



we obtain with probability one 



^2 x k el < oo, 



fc=0 



which together with the conclusion of Lemma A. 2 implies the almost sure uniform conver- 
gence of {X™} on Iq. 



FRACTIONAL BROWNIAN FIELDS OVER MANIFOLDS 



31 



□ 

As remarked in [T7] this result gives a sufficient condition for the existence of an almost 
surely continuous version of X x . The next result establishes Holder continuity. 

Theorem A. 4. (Manifold version of Thm 8.3.2 in [TJJ: Let the field X over Iq C M be as 
above and let 7 = sup{/3 : E[|Xr — X y \ 2 ] = o(d(x,y) 213 ) uniformly on Iq}. Then there exists 
a version of X with sample paths that are almost surely uniformly Holder continuous over 
Iq of any order (3 < 7. 

Proof. Let p = d(x,y). First note that, with p(u) as above, we have for any L > 

poo roo 

/ P(e~ x2 ) dx <c e / e -^-^ x2 dx < 00 

J L J L 

for any < e < 7 and some constant c e . But this is equivalent to 

diam(7 ) 

a/— log(n) dp(u) < 00. 

Thus by the previous result we have a version (which we also denote by X) for which 

\X X - X y \ < Bp(p) + C f P -\og{u) dp{u) a.s. 

Jo 

for some constant C > and some positive random variable B almost surely finite. 

Now for any < e < 7 we have some constant C t > such that p{p) < Cep 1-6 , and 
rp 

similarly / \J — log(it) dp(u) < C' e p' y ~ e for some C' e > 0. Thus, with probability 1, for 
^0 

each e > there is an almost surely finite positive random variable A e such that 

\X X -X y \ <A e d(x, y y- e Vx,ye/ . 

□ 

Note that we can also show under the hypotheses of the theorem that in any disk of 
positive radius in Iq the sample paths of X fail to be uniformly Holder of any order greater 
than 7. Indeed, 

X x — Xy 

d(x,y)-r+ e 

is a centered Gaussian random variable with variance 0(d(x,y)~z) and thus becomes 
almost surely unbounded as x — > y. For example we can pick any countable dense subset 
of Iq and modify A on a set of measure zero to obtain the failure of Holder continuity at 
each point in the set. Any stronger converse statement will require more refined tools, i.e., 
local times, which we will not attempt to develop here. 
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